The dielectric behavior of the living cell 
suspensions 

Camelia Prodan, Emil Prodan* 
University of Houston, 4800 Calhoun Rd, Houston TX 77204-5508 

08.04.1998 

Abstract 

In the limit of small concentrations and weak applied electric fields, 
the dielectric permittivity of suspensions of arbitrarily shaped, shelled and 
charged particles is calculated. It is proved that the dielectric behavior 
at low frequencies is dominated by the effects of the diffusion of the free 
charges on the shell surfaces. Our theoretical formula is valid in the 
low range of frequencies (a dispersion) as well as in the high range of 
frequencies (/3 dispersion). Will result that one can measure the membrane 
electrical potential by a simple investigation of the living cell suspension 
dielectric properties. 

1 Introduction 

The development of quantitative methods of characterizing the living matter 
represents the subject of many researches of our days. The study of the electric 
(dielectric) properties of biological systems belongs to this kind of programs. 
A simple measurement of the dielectric permittivity of a living cells suspension 
provides us with an important set of data which can be used to describe the 
living matter |^|, |l7j . There are many factors which may influence the dielectric 
behavior of the biological materials: structure, molecular dipoles orientation, 
self interactions, surface conductance, diffusion, membrane transport processes 
etc. All these factors influence one each other and it is hard (if not impossible) 
to separate the effect of a single one. However, some effects are dominant 
in certain ranges of frequencies and certain conditions. For example, in the 
/3-dispersion range (4 MHz-8 MHz), the dielectric behavior of the living cell 
suspensions is mainly influenced by the dielectric structure. At low frequencies, 
it is dominated by the diffusion effects of the free charges accumulated around 
the cell membrane J^]. Also, the orientation of molecular dipoles and membrane 
transport processes become relevant when strong electric fields are applied. For 
a survey of the field and an extended bibliography one can consult || . 
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The electric activity, in particular, the existence of the membrane potential 
is one of the simplest factor which allows us to discern between a living and a 
nonliving cell. We will show in this paper its influence on the dielectric proper- 
ties of biological materials. Our model do not include rotational effects or self 
interactions, so one expects that our results to be valid only for diluted suspen- 
sions and weak electric fields. We like to think that our work is complementary 
to the works done in electrophoresis, where the applied electric fields are very 
strong and the object of investigations is the induced membrane potential |p"lf . 
In this paper, the applied electric fields are weak and the object of investigation 
is the influence of the membrane potential. 

We calculate the dielectric permittivity of suspensions of shelled, arbitrar- 
ily shaped particles with a superficial distribution of free charges on the two 
faces of the shell, in the limit of small concentrations and weak applied electric 
fields. The method used here is similar with that of ||, the basic tools being 
the functional calculus and the spectral decomposition of the operators. For the 
real situation, the distributions of free charges are rather spatial than superfi- 
cial. For usual values of the (3-D) diffusion constants, they have the support in 
the immediate vicinity of the membrane (the distributions are practically zero 
for the distances larger than 1CP 9 m) and this allows us to treat them as su- 
perficial distributions. Our numerical application will reveal the dependence of 
the dielectric permittivity on the membrane electrical potential, mobility of the 
free distribution of charges and particle geometry. All these factors are equally 
important in the low range of frequency. It turns out that the consideration 
of non-spherical geometries (especially non-convex geometries) may play an im- 
portant role in some situations. For example, if one wants to investigate the 
living cell cycles by dielectric measurements on synchronized cell suspensions, 
then definitely, near the subdivision point, the shape effects are considerable [Q. 

We don't know yet how to relate our results with the electrorotation theory 
[ p~0| or how to combine them. We consider that is premature to compare our 
theoretical results with the experimental results obtained by electrorotation, 
which, at this time, seems to be the only accurate experimental data in the low 
range of frequencies. This is mainly because of different regimes of the applied 
electric field and concentrations. 



2 The general analysis 

From now, the perturbation of the external electric field (which will be called 
the excitation), due to the presence of a particle, will be called the response of 
that particle to the excitation. It is known that the Lorentz method of dielectric 
permittivity calculus works only if the response of the suspension constituents 
is linear in respect to the excitation. When diffusion effects are considered, this 
fact is no longer true. This is the reason we start our analysis with the response 
of a single particle to an external excitation. 

Let us consider a particle with two dielectric phases (shelled particle), placed 
in a spacial homogeneous and temporal oscillating electric field. In addition, we 
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consider that we have a superficial charge distribution on the shell faces. Laplace 
equation, satisfied by the electrical potential, has to be completed with the 
boundary conditions on the separation surfaces of the dielectric phases. These 
conditions will be deducted from charge conservation equation. If one considers 
a surface £ which separates two dielectric mediums, T>±, with a free charge p 
distributed on it, then the electrical current density is formed by a volume one, 

given by j ± = cr^ ■ E and a singular density, j , having the support on the 
interface £. Denoting with n the net charge concentration (note that a part of 
charges are due to different conductivities cr + , <J~), the continuity equation for 
a domain V, centered on the surface, is 



'^f ndv = ^ JtotaidS = j> j vo idS + j> i 

V dV dT> r 



dr, 



(1) 



where T = dV n £. Using the Maxwell equation, n = divD, it follows 



Jsjr. 



(2) 



dV 



The singular current is the only cause of the superficial distribution variation 



d_ 

dt 



P -dS= 4>j sg dr 



VnT, 



Finally, the passing equation of the electrical field, in the integral form, is 

d 



(a-E+-D).dS 



d_ 
~dt 



p dS, 



dV 



and in the differential form is 



n(a+ • E + + l£>+) n(a- ■ E~ + ^D~) 



dp 



(3) 



(4) 



(5) 



where n represents the normal at the interface. For shelled particles, the com- 
plete system of equations will be 



A$ = ; x e 3fJ 3 \(Ei U S 2 ) 

-»+ d -+ -. 
n(a+ -E +-^D ) - n(a~ 

x e Ei U E 2 

div s j sg = -^7; ieSiUE 2 



dt ' dt 



(6) 



E 



dt 

Eq ■ cxp(ju>ot) as 



where £1, £2 are the external respective the internal face of the shell and 
E exp (ju>ot) is the external field. 
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The temporal Fourier decomposition leads to the following form of the equa- 
tions and boundary conditions: 



A$ = 0;ie 5R 3 \(Si UE 2 ) 
(cr + + ju;e + )n ■ E + - (cr~ + ju;e~)n ■ E~ 
x e Ei US 2 

div v j sg = -jup; x e si u s 2 
J £ ; w = w 

! 0; w ^ w 



(7) 



as 



CO. 



For w^O, the equation e?iv E j = —jup imposes 
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p(u)dS = ~ J div s j sg dS = 0. 

Si 2 



(8) 



This means that for ui cj and w 7^ 0, we obtain trivial boundary conditions. 
So the system has a trivial solution. For w = 0, the above equation cannot be 
written. It is being replaced with the condition 



p(0)dS = ±Q, 



(9) 



S l,: 



where ±Q represent the total free charge distributed on the two interfaces. 
This leads to nontrivial boundary conditions and consequently to a nontrivial 
solution. In consequence, the solution of the system has the following general 
form: 



®(x,t) = $ (x) + ${x) ■ exp(jw i) 
p(x, t) = p Q (x) + p(x) ■ exp(jw i) 



(10) 



The components $0 and p represent the solution of the system for zero external 
excitation and will represent the equilibrium of the system. The equilibrium will 
be discussed in the next section. The components $ and p are solutions of the 
following system: 



A$ = ; x e 5R 3 \(Si U £ 2 ) 

(a+ +juje+)-— - (a +juje )— - 
_^ On On 

div s 3 sg = -jup; x e Si u £ 2 

E — ► E as : - - x . 



: jup; ieEiUS 2 



(11) 



where the pulsation of the excitation was re-denoted with u>. Introducing the 
expression of the singular current 



j sg = J conduction + j diffusion = ~7i 1 V Ei $ - AVs.ft, i = 1,2, (12) 
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and working with the complex electric permittivity e* = e + a/ (jco), we can 
write 



' A$ = 0;i£ K 3 \( s i U S 2 ) 

on _ an _ (13) 

E — > Eo as x — > oo. 

where the index £ means that the operators are calculated on the surface. The 
system is not linear, because the conductivities of the superficial charges, -f i , 
depend on p. We will consider in the following that we are in the limit of weak 
external electric fields so we can consider that the conductivities of the free 
charges are given by the equilibrium configuration. Thus the system becomes 
linear. Now, the total electric field will be the sum of the electrical fields of 

the equilibrium configuration and of the perturbation, E to tai — E e + E. In 
the weak fields approximation, the average on the different orientations of the 
particle leads to the following value of the electric field inside of the particle: 

E) = E - — [e~ -N d£l~, (14) 

/ 47T / N AT V ; 

/ orientations J 



where E-^ is the electric field (only the perturbation part) inside of the particle, 

when the external electric field, N, is of norm one. Further, this average can be 
calculate with the formula 

e) =±(V\v_;r x )£„. ,ir» 

/ orientations 

where i designates three orthogonal directions. The average electric field is 
proportional with the excitation. This allows us to use the Lorcnz formula 

£ sus = So (l + 1 (16) 

for the suspension dielectric permittivity calculus. Here p is the volume concen- 
tration of the suspension and a is the polarization of the particles 

a = ^-[dv £ —^(E) (17) 

y * ' orientations 

where V is the particle volume, e c is the dielectric permittivity of the exterior 
medium and e is the dielectric permittivity of the particle. 
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3 The equilibrium 



In the following the primitivities and the conductivities will be denoted by: so, 
Co for the exterior, e\, <j\ for the membrane, £2, cr 2 for the interior. The sys- 
tem which establish the equilibrium configuration is obtained from the general 
system by cancelling the temporal derivatives: 

A$ e = ; x e 5R 3 \(S! US 2 ) 

-=- bi= °i— -=- ; 1 = M , . 

dn <9n (18) 

J sg = 

Je 1|2 ^ = ±Q. 

If we consider the conductivity of the membrane, <7 1; equal with zero, we obtain 
trivial Neumann condition for the electric potential, inside and outside of the 
particle. This imposes constant electrical potential in the two regions. The 
distribution of charges on the two interfaces, n^, can be calculated by applying 
the Gauss law: 

T i = £i— — W, T2 = £1—^ |s 2 (19) 
on on 

where r includes both, the free charge p and the charges which are accumulated 
on the interfaces due to the conduction currents. The equations for free charges 
are 

divfrt ■ E t - £>Vp J = -DAp 0i = p . = const., (20) 

where £? t is the tangent electric field to £1,2- The link between p and Q results 
by the reason that the inside region is not an electric charges reservoir, which 
means that the total charge on the interfaces is given only by the free charges. 
Thus 

±Q = J P 0l JS = p 0l2 ■ S h2 ; Q = C ■ AV (21) 

Sl,2 

where AVo is the electric potential difference between the two surfaces (mem- 
brane potential) and C is the capacity of a condenser having the geometry of 
the shell and dielectric permittivity equal with z\. The values of the p . fix the 
electric conductivity on the two interfaces 

C 

7i,2 =Po 12 -"1,2 = ±"1,2-5— " AV °> ( 22 ) 

<->l,2 

where Uj, i = 1, 2 are the charge mobilities on the two interfaces. 

It was seen that, in the calculus of the dielectric permittivity, only the per- 
turbation part of the electric field is important. So the only influence of the 
membrane electric potential on the dispersions curves of the dielectric permit- 
tivity comes from the above formula. 
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4 The effective calculus of the polarization 



The charges on the second interface lie behind the shell, so the external electric 
field will have a smaller influence on them than on the charges of the first 
interface. In the same time, it is an experimental fact that the mobility of 
the inside charges is much smaller than that of external charges. Both reasons 
justify the idea of considering the free charges of the inside interface to be 
fixed (i.e. zero mobility). With this simplification, the value of p 2 is fixed at 
the equilibrium value. Now we have inhomogeneous particles with free electric 
charges only on their surface. 

4.1 The equivalence between an inhomogeneous particle 
and a homogeneous one 

Let us consider first the situation of a homogeneous particle, with dielectric 

permittivity e and conductivity a, placed in the electric field E^e^ 1 . If we use 
the expression of the simple layer for the electric potential 

$(x) = -x-E + -L ■ / ^ V l dS y , (23) 



47T 

yes 



x - y 



the passing equation for the electric field through the surface of the particle, E, 
becomes the following integral equation for the charge distribution /x: 

^(x)-E[ t j](x) = n-E , (24) 

e* - e* 

where A = and E is the operator 

e*+e* Q 

E M (*) = ^ • J {x - y) J x n{y)dS y . (25) 



yes 



x - y 



Using the spectral decomposition of the E operator, the solution of this equation 
is 

M = Et— T Pn[n-E ], (26) 

„ 2 ^' Xn 

where P n is the spectral projector corresponding to the eigenvalue x n - The E 
operator is not a symmetric one, so some precautions are needed. About this 
operator and the validity of the above decomposition, one can consult p5| . 

Let us consider now the shelled particle placed in the same electric field. We 
consider here that the external surface of the shell is obtained by expanding 
the internal surfaces by a factor 5 > 1. The equivalence problem between this 
inhomogeneous dielectric particle and a homogeneous one can be formulated in 
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the following terms: there exists a distribution of charges only on the external 
surface of the shell, fi e , which provides the true electric field outside of par- 
ticle. In the dipole approximation, this distribution is given by the following 
expression: 

ft=El-^ Pn[n-E ]. (27) 

where A„ are given by: 

\ _ £ n ~ £ o _* _ _* X" + 2 + ^ ( 2Aj ~~ Xn) . _ ^2 ~ £ 1 /r,n\ 

Comparing with the formula for the homogeneous particle, we conclude that 
the shelled particle is equivalent to a homogeneous particle which has the same 
geometry and volume. The equivalation relation is given by 

£ * - {<}» , or e* = $>*P„ ~* E e ^™- ( 29 ) 

n n 

In the particular case of the sphere, the above formula is identical to the existent 
formula jl2| . We will give in the following a sketch of the proof. For the shelled 
particle, the single layer expression of the potential is 

$ (a:) = -xE + 



1 

4-7T 



My) 



x - y 



dSy 



1 

47T 



yes 2 



x - y 



(ISy 



(30) 



The electric potential created by /x e will be denoted with <f> e . The potential <E> e 
is equal outside of particle with the exact potential, $, only if the two potentials 

fulfill the same Neuman boundary conditions: ^ = . Using the single 
layer expression, the above condition can be written as 



~ E W] = \v\ - E u M - E 12 [M 2 ] 



where 



47T 



(x - y) ■ n x ^ 



yes j 



x - y 



(31) 



(32) 



The passing equations of the exact electric field, $, through the shell surfaces 
take the form 



2x7^ - Eu [fi{\ - E 12 [n 2 ] = nEo 
aijM2 - E 2i [Mi] ~ #22 [M2] = 



(33) 
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where = -j j— Here comes the dipole approximation: 

#12 M ~ ^#22 D*a] (34^ 

which says that the normal component of the electric field created by the charge 
distribution p 2 on the surfaces Si has the same angular dependence as on the 
surfaces E 2 but the strength of the field is diminished by a factor S~ 3 , and that 
the normal component of the electric field created by the distribution p 1 on the 
surface £2 is equal with the normal component on the surface Si. With the 
observation that the operator E is scale invariant, the above equations become 

I 2lT^i _ E K] - jgE \p 2 ] = nE ^ 
\ 2X^2 - E [/ij - E \p 2 ] =nE 

Using again the spectral decomposition of the operator E and the link (31) be- 
tween p e and /j,i , /i 2 , the desired form of the distribution p e follows immediately. 
One can see that the method can be extended to the multi-shelled particles. In 
this way, the results of Q can be extended to arbitrary geometries. 

Now, we want to consider the situation when the free distribution of charges 
p is present. For a homogeneous dielectric particle with dielectric permittivity 
£ and conductivity a which has a free charge on its surface, the equations are 



• E + \p + E[fj]J - e* ( n ■ E Q - \n + E[fi] ) = p 



(36) 



where As is the Laplace operator on the particle surface. Here we have contin- 
ued the idea of the second section by introducing the simple layer expression of 
the electric potential. We can pass to the situation of a shelled particle by using 
the equivalence relation e* — > ^2 n e^Pn- (Here we use our assumption that the 
distribution of the free charges of the internal shell surface can be considered 
immobile). The equations become 

n ■ E Q + \p + E[p] 
~ (E< • A>) (n>Eo- \p + E\p] \ = p ^ 

(A S r + ^)p=-*A Sl $ 

Here, p denotes the free distribution of charges of the external shell surface. 
One observes that the single layer distribution can be decomposed into p = 
Pp + p a , where p^ solves the equations for p = 0. It turns out that p^ is 
practically constant at low frequencies and that p a is practically zero at high 
frequencies. In consequence, the pg distribution is responsible for the behavior 
at high frequencies (the range of the beta dispersion) and the addition of p a 
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will give us the behavior at low frequencies (the range of alpha dispersion) . The 
equations for jj, a and p will be: 

[ E (k( £ n + £ o) - Xn(£n ~ £ *o)) Pn H a = P 

i n loo) 

where ( $ a ) represents the electric potential generate by the p^ (p, a ) dis- 
tribution. The spectral decomposition of the Laplace operator leads to the 
following equation for \i a : 

E (U £ n + £ *o) - XnK - £*o)) ^ P a 



» Si 



where Gi is the spectral projector of the Laplace operator corresponding to the 
eigenvalue Denoting with if the operator 



(41) 



the equation takes the form: 

{fE(i« + ^)-x„«-<))4 

We will discuss later our method of solving this equation. 
4.2 The formula for polarization 

We can use the idea of equivalence to find a compact formula for polarization. 
For this, we have to transform the volume integral which appears in the defini- 
tion of a, 

a = 1 f dSl^fdv £ —^N ■ E, (42) 



V J N 

V 

into an integral on the external side of the particle surface where we can use the 

equivalence. In the above formula, e depends on the point, and E is the exact 
field inside of the particle, prodeced by the excitation N, of norm one. We will 
calculate the term as follow: 

/ dv ■ -E = Jdv ^grad • ivj • 3 + div j 

= J_ Jdv div (^(x ■N)-D S j+ Jdv div (n$J (43) 

= — J (x ■ N)D ■dS+ J$N-dS 
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where, if it is used the passing relations for D, D+ — D n = p and the continuity of 
the potential, the integrals can be processed on the external side of the surface. 



fdv £ —^N.E=±f(x-N)(-p + D+)-dS 

+ / <S>N • dS = — J(x-N)(-p + e Q E+) -dS+ f$N- dS. 

S+ £ °E+ E+ 



Outside of particle, the electric field is equal to that of the equivalent homo- 
geneous particle. The normal component of the electric field can be expressed 
with the help of the E operator: 

d$ + — — 1 
En = — -=r = n ■ N + -fj, + Ep, (45) 
on L 

where /x is the equivalating charge distribution (solution of (37)), resulting: 

j- dv ■ E —^N -E = — J(x-N) (-p + e (n-N+±p + Ep\) dS 
-JN .E e dv=-f(x-N) ( £e* ■ P n )(n ■ N - (± - Xn )p) dS (46) 

V £ o s - \ n J \ J 

- J (x ■ N)E e dS 



E- 



Having that nE e = ~7pr = n ■ N — ^p + E [p], we can continue 



J dv . S -^N ■ E 

V e ° 

= /(* ■ N) (e^^ • A.) (n ■ N - (I - Xj >) dS (47) 
= E^-^ • Rx ■ N)P n (n-N-{\- Xn )p) dS. 

If one denotes the scalar product on the L 2 (E) functions space by ( , ), 

(<P, <!>) = j ^(x)<P(x) dS Xl V if, <j> e i 2 (S), (48) 

E 

we obtain 



Jdv ■ - — —N ■ E = x 



n <^o 



^x-N\P n \n-N^-^- Xn )(^x-N\P n \p^ 



(49) 
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Having in view the discussion about the average on the orientations, the angular 
integral can be eliminated and the final formula for polarization becomes 



3V< 



« = wE 1 ^ \(x-Ni\P n \n-Ni 



I _ \1 (50) 

-{\-Xn)(x-Ni\P n \ljj 

where i = 1,2,3 denotes three orthogonal directions and [i % denotes the total 
distribution of the single layer corresponding at the excitation JVj. If the free 
charge is not present, the above formula becomes 

" = 3fE (£ "~ £;) {(x-Ni \P n \n-Ni 

n,i I \ 

-(l-xJ^-^lTZ^A.In-TV^J (51) 
= wE i _ A " r (x -Ni\ P n \n -Ni 

This formula is exact in the high frequencies range and it describes correctly 
the f3 dispersion region. The exact formula of a must be completed with the 
term which comes from fi a 

a - w£ i _ A " i (x-Ni\ P n | n ■ N y 

n,i 2 Xn^n \ . r ., . 

(T _ y V e * _ e *) / „ \ ^ 0Z J 

n,i \ 

As we already discussed, the last term is practically zero at high frequencies 
and becomes dominant at low frequencies. 



5 Numerical results 

First, we will present our method of solving the complicated equation (41) of 
fi a . Here, the principal problem is the calculation of the spectral projectors of 
E and As operators. Let us suppose that our surface is given in the spherical 
coordinates by r — r (9, (p). To solve the problem for the E operator, we will 
choose the following orthonormal basis in L 2 (£): 

OW mez ; y lm (0, <p) = (53) 
l>\m\ ^/<r(6,ip) 

where {Yi m } l m represents the spherical harmonics and dS = a (9, ip) sin (9) d9d(p. 

The first step is the calculation of the matrix elements of the operator E. Then, 
taking acceptable truncated matrices, we will calculate the approximative spec- 
trum. The approximative projectors will be given by P n = P n (^E^J , where P n 
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are the polynomials which satisfy P n (x m ) = S nm for each eigenvalue x m - For 
the operator As, we have followed the way presented in Q. 

We have chosen the surfaces of figure 1 for our numerical application. As 
we already emphasized, 7 is proportional with the electrical potential of the 
membrane and the mobility of the free charges. Because the diffusion constant 
is also proportional with the superficial charge mobility, the quantity j/D is 
proportional only to the membrane electrical potential. We will choose different 
values for the diffusion coefficient, D, and for the ratio j/D, which will give 
us the dependance of the dispersion curves on the mobility of free charges and 
membrane electrical potential. Because the a effect on the conductivity is much 
smaller than on the dielectric constant, we are enforced to consider two different 
values of the parameters: one for which one can see in detail the behavior 
of dielectric constant and one for which one can see in detail the behavior of 
conductivity. The results are presented in figures 2-13. The following parameters 
were kept constant in our numerical analysis to the following values (the units 
are in S.I): e = 78e mc , £1 = 10e„ QC , e 2 = 50e vac , &o = 0.2, o\ — 0, er 2 = 0.3, 
p = 0.06, 5 = 1.004. 



6 Conclusions 

We want to point out first the relation of our results with the existent results. 
In the /3-dispersion range of frequencies, the importance of E operator, in this 
type of analysis, was pointed out in p8[ . Basically, this method reduces the 
original equations of the electric field to surface equations. Mathematically, this 
means that it transforms a problem involving unbounded operators (Laplace 
operator corresponding to certain boundary conditions) to a problem involving 
compact operators (the operator E). Here and in Q it is given a semi-empirical 
justification of the formulas (28) but, unfortunately, the formula of polarization 
is not exact. Our paper states in very precise terms what equivalation means, 
proves rigorously the formulas (28) and that of polarization. Another result 
directly related with our approach would be , where the method of equivala- 
tion is discussed. Note that this result (and the method itself) depends on very 
particular surfaces considered there. In the a-dispersion range of frequencies, 
there are some results for the case of spherical geometry || (and what follows 
after) but, as is stated in the beginnig of this paper, the results are valid only 
in a quasi-statical regime. One can see that, in the limit of zero diffusion con- 
stants, its passing equations does not reduce to the usual equation, but rather 
to two separate equations: one for conduction current and one for displacement 
current. In fact, this was our reason for we have treated very carefully the 
introductory sections. 

In conclusion, this approach gives a unitary treatment of a and /3-dispcrsion 
in the general context of arbitrary geometries. In the /3-dispersion range of fre- 
quency, our paper corrects and puts some existent results in rigorous settings. 
In the a range, it proposes a completely new approach which allows us to avoid 
problems like that mentioned before. Also, it proves the general formula of po- 
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larization which includes the effects of the free charges. Our model reproduces, 
qualitatively, the dielectric behavior of living cell suspensions in both a and (3 
frequencies ranges and is flexible enough to reproduce any given curve which 
has this shape. The numerical application shows a very strong dependency of 
the dielectric dispersion curves on the membrane electrical potential. This fact 
anticipates the possibility of the experimental measurement of the membrane 
electrical potential by a simple measurement of the dielectric permittivity. Be- 
sides this important cell parameter, one can measure the mobility of the free 
charges and, as in |T^| , the volume concentration of the suspension. We em- 
phasize that this method is a non-destructive and can be a very fast one (a 
discussion about this can be found in ||). 
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Figure captions 



Fig. 1 The surfaces chosen for our numerical application: 
r{9) = (l + acos0)/(l + a), where a = 1,2,3. 

Fig. 2 The dcpendance of dielectric permittivity on membrane potential. 
The numerical values are a = 1, D = 10~ 8 , 7/D = 1 (plus), 7/D = 3 (circle) 
and 7/D = 5 (star). 

Fig. 3 The dependance of conductivity on membrane potential. The nu- 
merical values are a = 1, D = 1CP 8 , 7/D = 1 (plus), 7/D = 3 (circle) and 
7/D = 5 (star). 

Fig. 4 The dependance of dielectric permittivity on membrane potential. 
The numerical values are a = 1, D = 1(T 8 , 7/D = 0.1 (plus), 7/D = 0.5 (circle) 
and 7/D = 1 (star). 

Fig. 5 The dependance of conductivity on membrane potential. The nu- 
merical values are a = 1, D = 10~ 8 , 7/-D = 0.1 (plus), 7/-D = 0.5 (circle) and 
7/D = 1 (star). 

Fig. 6 The dependance of dielectric permittivity on mobility. The numerical 
values are a = 1, 7/D = 5, D = 10~ 7 (plus), D = 10~ 8 (circle) and D = 10~ 9 
(star). 

Fig. 7 The dependance of conductivity on mobility. The numerical values 
are a = 1, 7/D = 5, D = 10~ 7 (plus), D = 10~ 8 (circle) and D = 10~ 9 (star). 

Fig. 8 The dependance of dielectric permittivity on mobility. The numerical 
values are a = 1, 7/D = 0.5, D = 10~ 7 (plus), D = 10~ 8 (circle) and D = 10~ 9 
(star). 

Fig. 9 The dependance of conductivity on mobility. The numerical values 
are a = 1, 7/D = 0.5, D = 10~ 7 (plus), D = 10~ 8 (circle) and D = 10~ 9 (star). 

Fig. 10 The dependance of dielectric permittivity on geometry. The nu- 
merical values are 7/D = 0.5, D = 10~ 8 , a = 1 (plus), a = 2 (circle) and a = 3 
(star). 

Fig. 11 The dependance of conductivity on geometry. The numerical values 
are 7/D = 0.5, D = 10~ 8 , a = 1 (plus), a = 2 (circle) and a = 3 (star). 

Fig. 12 The dependance of dielectric permittivity on geometry. The nu- 
merical values are 7/D = 5, D = 10~ 8 , a = 1 (plus), a = 2 (circle) and a = 3 
(star). 

Fig. 13 The dependance of conductivity on geometry. The numerical values 
are 7/D = 5, D = 10~ 8 , a = 1 (plus), a — 2 (circle) and a = 3 (star). 
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